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Abstract 

The Standard Model provides an excellent description of the observables measured at 
high energy lepton and hadron colliders. However, measurements of the forward-backward 
asymmetry of the bottom quark at LEP suggest that the effective coupling of the right- 
handed bottom quark to the neutral weak gauge boson is significantly different from the 
value predicted by the Standard Model. Such a large discrepancy may be the result of a 
mixing of the bottom quark with heavy mirror fermions with masses of the order of the 
weak scale. To be consistent with the precision electroweak data, the minimal extension 
of the Standard Model requires the presence of vector-like pairs of SU{2) doublet and 
singlet quarks. In this article, we show that such an extension of the Standard Model is 
consistent with the unification of gauge couplings and leads to a very rich phenomenology 
at the Tevatron, the B-factories and the LHC. In particular, if the Higgs boson mass lies 
in the range 120 GeV < rrih ~ 180 GeV, we show that Run II of the Tevatron collider with 
4-8 fb^^ of integrated luminosity will have the potential to discover the heavy quarks, 
while observing a 3-0" evidence of the Higgs boson in most of the parameter space. 



1 Introduction 



In the absence of direct evidence for physics beyond the Standard Model (SM), precision elec- 
troweak tests are the best way to get information about the scale and nature of a possible 
breakdown of the SM description. While the SM has held firm in the face of a great number of 
precision electroweak tests, the model has not emerged completely unscathed. Fits of the SM 
to electroweak data show about a 2.5-0" deviation in the b-quark forward-backward asymmetry 
(A^p^) and this situation has not improved much in the last five years. This discrepancy is 
important for two reasons. On one hand, it seems to indicate a significant deviation of the cou- 
pling of the right-handed bottom quark to the Z-gauge boson (see, for example, Ref. 0). On 
the other, this measurement plays an important role in the present fits to the SM Higgs mass; 
the removal of the heavy quark data from the electroweak fits would push the central values of 
the Higgs mass to lower values, further inside the region excluded by the LEP2 searches 

There are two ways of solving this apparent discrepancy, and both of them seem to indicate 
the presence of new physics. In Ref. ^ it was proposed to exclude the heavy quark data while 
introducing new physics that raises the central value of the Higgs boson mass, and improves the 
fit to the other observables. Such a task requires new physics that gives a negative contribution 
to the S parameter, positive contributions to the U parameter and a moderate contribution 
to the T parameter. At least two examples of this kind of physics have been presented in the 
literature |3], [S]; the first within low energy supersymmetry and the second within a warped 
extra-dimension scenario. 

An alternative to this procedure is to take seriously the heavy quark data while introducing 
new physics that modifies in a significant way the right-handed bottom quark coupling to the Z. 
The Beautiful Mirror model of Ref. |H] accomplishes this by allowing the b-quark to mix strongly 
with a set of exotic vector-like quarks. This model turns out to have several other interesting 
features which we investigate in this paper. To be specific, we consider the unification of gauge 
couplings, additional patterns of flavour mixing, the Higgs phenomenology, and searches for 
the heavy vector quarks. 

The model consists of the SM plus additional vector-like "mirror" quarks. These are a 
pair of SU{2) doublets, \I'l,_r = C^f'^), and a pair of SU{2) singlets, ^. Here and in what 
follows we use primed fields to denote gauge eigenstates, while mass eigenstates are written as 
unprimed fields. The gauge group quantum numbers are the same as those of the analogous 
SM particles: (3, 2, 1/6) for the doublets, and (3, 1, —1/3) for the singlets. Since the quarks are 
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added in vector-like pairs, these can have gauge-invariant Dirac masses, and the model is free 
of anomalies. This is a minimal set of mirror quarks needed to improve the fit to electroweak 
data. 

The Yukawa and mass couplings of the mirror quarks are taken to be 

- (2/3QL + yA)i'R^ - M^Wr + {h.c.) 

where Q'^ = {^y) is the usual third generation SM quark doublet, and $ = (^0) is the Higgs 
doublet. This is the most general set of renormalizable couplings provided the mirror quarks 
couple only to each other and to the third SM generation.^ As pointed out in [6 , the Yukawa 
couplings yh: ys and y^ are constrained to be much smaller than ?/2- Adjusting the ratio 
{-^y2)lMi ^ 0.7, where v = 246.22 GeV is the Higgs VEV, gives the best fit to precision 
electroweak data while reducing the discrepancy in to about one standard deviation, and 
keeping the left-right b-quark asymmetry within one standard deviation of the value measured 
at SLC This forces ?/2 to be 0{1) since Mi > 200 GeV is needed to explain why mirror 
quarks have not yet been observed. On the other hand, there are no strong constraints on y^, 
which mixes the exotics and therefore has only a small effect on the SM sector of the model. 
Following jH], we will mostly neglect y^ for simplicity, altough we will comment on the effects 
of this coupling whenever they are relevant. 

This paper consists of seven sections. In Section 2 we examine the running of the gauge 
couplings and their unification at a high scale. In section 3 we discuss the issue of flavour 
mixing as well as the quark couplings to the neutral and charged weak gauge bosons, and the 
Higgs. Section 4 consists of an investigation of the Higgs phenomenology in the model. In 
Section 5 we review the current limits on exotic quarks and investigate the possibility of finding 
mirror quarks at the Tevatron. In Section 6 we examine how the new types of flavour mixing 
possible with mirror quarks can affect CP violation in B —>■ (pKg decays. Finally, Section 7 is 
reserved for our conclusions. 

2 Unification of Gauge Couplings 

The idea that the low energy gauge forces proceed from a single grand unified description is a 
very attractive one, and is supported by the apparent convergence of the weak, hypercharge and 
^Note that couplings like Q'^^'j^ and can be rotated away. 
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strong couplings at short distances. The interest in low energy supersymmetry, for instance, 
has been greatly enhanced by the discovery that the value of the strong coupling, as{Mz), 
can be deduced if one assumes that the gauge couplings unify at a high scale. This prediction 
depends on model-dependent threshold corrections at the GUT scale, but to within the natural 
uncertainty in these corrections JT], the predicted value of as(M^) is perfectly consistent with 
the values measured at low energies. In the Standard Model, instead, the assumption of gauge 
coupling unification leads to a prediction for as{Mz) that differs from the measured value by 
an amount that is well beyond the natural uncertainties induced by threshold corrections. 

In jHj it was noted that, to one- loop order, adding mirror quarks of the type considered 
here to the SM greatly improves the prediction of as{Mz) based on the assumption of gauge 
coupling unification. We extend this analysis by including the two-loop contributions to the 
gauge coupling beta functions and the low-scale threshold corrections. Since, for the consistency 
of this study, the Higgs sector must remain weakly-coupled while the Higgs potential should 
remain stable up to scales of the order of the unification scale, Mq, we also investigate the 
related issues of stability and perturbative consistency of the Higgs sector. 

In extrapolating the low energy description of the theory to short distances, it is important 
to remark that the Beautiful Mirror model [0] does not provide a solution to the hierarchy 
problem. Therefore, a main assumption behind this extrapolation is that the physics that leads 
to an explanation of the hierarchy problem does not affect the connection of the low energy 
couplings to the fundamental ones. An example of such a theory construction is provided by 
warped extra dimensions jH], and has been investigated by several authors [ni-[II]- In order 
to preserve the good agreement with the precision electroweak data, the Kaluza-Klein modes 
must be heavier than a few TeV in this case ^2]; and therefore the low energy physics analyzed 
in the subsequent sections will not be affected. On the other hand, extra dimensions could 
modify the proton decay rate in a significant way by introducing new baryon number violating 
operators, and, in the case of warped extra dimensions with a Higgs field located in the infrared 
brane, would make the issue of the running of the Higgs quartic coupling an irrelevant one. For 
the rest of this section, we shall proceed with a pure four dimensional analysis of the evolution 
of couplings and of the proton decay rate. 
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2.1 Renormalization Group Equations 



Using the results of OCl], the two-loop {MS scheme) gauge coupling beta functions are 

where t = In (y^^ is the energy scale, and / = 1,2,3 refers to the U(l), SU(2), and SU(3) 
gauge groups respectively. The first term is the one-loop contribution, while the other terms 
come from two-loop corrections. 

The coefficients hi and hki are given by 

6i = -i, 62 = 1, 63 = 5, (3) 

and 
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(4) 



In the SM, the corresponding one-loop beta function coefficients are hf^^ = —41/10, fof*''^ = 19/6 
and 6f*^ = 7. The variation of these coefficients are hence Abi = 2/5 and A62 = A63 = 2. 
Since 62 and 63 are shifted by an equal amount, they tend to unify at the same scale as in the 
SM, about a few times 10^^ GeV. Interestingly enough, the shift in bi is much smaller than 
that of 62 and 63, leading to, as we shall see, a successful unification of the three couplings. 

The coefficients Yj^F) involve the Yukawa couplings. Neglecting the small Yukawa couplings 
yb,y3,y4, and y5, they are 

Yl{F) = Cityt' + Ci2y2' (5) 



where 
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Cif=\ i 11 (6) 
2 2 



and f = t,2. 

The Yukawa couplings evolve according to 

where f = t,2. The one-loop and leading two-loop contributions to jSf were calculated follow- 
ing Of the two-loop terms, we include only those involving or the Higgs self coupling 
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A; the terms are enhanced by large colour factors while the A terms can become important 
when investigating the stability of this coupling. The one-loop contributions are 

(4vrM^ = IvI + Sy^-i^lgf + lgl + Sgiy (8) 
The two-loop contributions that we have included are 

(4vrr/?f = h'-6ylX + gl{20y^ + ^6yl)-^gl (9) 

The total beta-function is the sum of these pieces: (3/ = j3^p +(3j'\ Aside from the modifications 
due to the mirror quarks, these are in agreement with the results of [T^ . 
The Higgs self-coupling A is taken to be 

^ D /i2<l>t$ - iA($t$)^ (10) 

With this definition, the tree-level Higgs mass is vrih = y/Xv, where v = 246.22 GeV = 
-\/2 (0°). A evolves according to ^ = Px. We have calculated the one-loop and leading two-loop 
contributions to Px using the results of [13 E] . As for the Yukawas, only the largest two-loop 
terms involving or A were included. For the one- loop part, we obtain 

(4vr)^/5« = 12A^ - 0,? + 9,,^) A + ^ + ^glgl + g^^ 

+ 12X{y^ + yl)-12{yt + y'^). (11) 

The two loop part is given by 

(47r)^/3f = -78X^-72iy^ + yi)\'-3iyt + yl)X + 60iy^ + yl) 

+ 18{^gl + 3gl)X^ + 80gl{y^ + yl)X - QAgHyf + y',). (12) 

Again, the total beta-function is the sum of the one- and two-loop parts. 
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2.2 Input Parameters and Threshold Corrections 



We have investigated the running of these couphngs numerically. The initial values MS scheme 
values were taken from [T^ : 

a-^{Mz) = 127.922 ±0.027 sin^ 9^ (Mz) = 0.23113 ± 0.00015 

Mz = 91.1876 ±0.0021 v = 246.22 GeV (13) 

rhtimt) = 165 ±5 GeV. 

These parameters correspond to the effective SU{3)c x U{l)em theory with five quarks obtained 
by integrating out the heavy gauge bosons and quarks in the full SU{3)c x SU{2) x U{1)y theory 
at scale Mz- Threshold corrections to the gauge couplings arise in the process of matching the 
theories. We define 

«r' = |(l-sin2^^)«-\ (14) 
5 

6l2^ = sin^ 6'^q;^\ 

Then the gauge couplings at Mz are given by 

«z~^ = «r^ + Pi> (15) 
where the pi terms represent threshold corrections. To one-loop order, these are [7j 

P3 



^yin(^ 

^M, 

C 



P2 = sin^ 9^ 



bn Mz vr ^ ^ Mz 



3 /I - sin^6'^\ 

where the sums run over ( = cj,^. As shown in j6j, the tree-level masses of the mirror 
quarks are given by 

= Ml, = iMf + Yif\ 

= M2, ^ ^ 

where Y2 = ■^y2- These parameters are not completely independent. As explained above, 
Y2 — 0.7 Ml [6j gives the best fit to electroweak data, while M2, Mi > 200 GeV are needed to 
explain why these exotics have not yet been observed at the Tevatron [TH] (see section 5). 
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Figure 1: Range of as consistent with a 1% unification of tlie gauge couplings plotted against 
the unification scale Mq. 

2.3 Numerical Evolution 

The unification of gauge couplings was investigated by fixing sin^ 6y^{Mz) and aem(^z) accord- 
ing to their measured values, and varying as{Mz) until the gauge couplings unified to within 
1%. GUT-scale threshold corrections were not considered. Figure [T] shows the range of as{Mz) 
obtained in this way for 250 GeV < M2 < 1000 GeV and all values of \{Mz) and 1/2 (M^) 
consistent with unification. (See the following section.) The range is plotted against the uni- 
fication scale. In general, the unification is quite insensitive to the input values of M2, A, and 
?/2- The scale of unification is quite high, Mq = (2.80 ± 0.15) x 10^^ GeV, depending on the 
input values, at which point the unified gauge coupling constant has value = 35.11 ± 0.05. 

The predicted range of the strong gauge coupling is in excellent agreement with the values 
measured experimentally. This agreement is quite intriguing since the particular completion 
of the Standard Model considered in this work is motivated by data and not by any model 
building consideration. We shall not attempt to construct a grand unified model leading to 
the appearance of the mirror quarks considered here in the low energy spectrum. Instead, 
we will concentrate on additional issues regarding the renormalization group evolution of the 
dimensionless couplings of the theory, as well as on exploring general features of the low energy 
phenomenology of this particular extension of the Standard Model. 
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2.4 Stability and Non- Triviality of the Higgs 

If the extrapolation of the model up to high scales is to be self-consistent, it should remain stable 
and weakly-coupled up to the unification scale. The only source of trouble in this regard is the 
Higgs self-coupling A. Stability of the Higgs sector requires X{Q) > 0, for all Q < Mq while 
perturbative consistency means that A must not be too large. For concreteness we demand ^ 
that < A < 2 up to 10^'' GeV. This is sufficient to guarantee that the effective potential is 
similarly well behaved [20] • The evolution of A is largely determined by the initial values of ?/2 
and A. Only for a small subset of initial values does A remain well-behaved (i.e. < A < 2) up 
to Mq- This subset is shown in Figure El where we have written A in terms of the (tree-level) 
Higgs mass, and Y2 in terms of assuming Y2 = 0.7 Mi [Hj. 

We compare the allowed region with the region favoured by precision electroweak data found 
in . There is a small overlap between the allowed band found here and the 1-a allowed region 
of corresponding to 160 GeV < < 180 GeV and < 225 GeV (m<^ < 275 GeV). We 
show in sections 4 and 5 that these mass ranges will be probed by Run II at the Tevatron with 
4-8 fb"^ of data. 

Finally, we note that including the coupling (see Eq. ((T))) would tend to displace the 
shaded region in Figure El downwards, slightly increasing the preferred range of Higgs masses. 
This is because, to one loop order, the beta functions are modified as 1/2 — > (y| + y|) by this 
inclusion, while the best fit value of Y2/M1 changes very little (see section IXT|) . Thus, our 
bound on 1/2 obtained with = translates into a bound on \/y\^ y\ in the more general 
case, leading to lower values of for a given rrih- 

2.5 Proton Decay 

Grand unified models induce baryon number violating operators that lead to a proton decay 

rate that may be observable in the next generation of proton decay experiments. The present 

bounds on the proton lifetime [21] already put relevant bounds on grand unified scenarios. In 

four dimensional supersymmetric grand unified models, for instance, dimension five operators 

may easily induce a proton lifetime shorter than the present experimental bound This 

situation may be avoided by a suitable choice of the low energy spectrum pHj. Heavy first and 

second generation sfermions and light gauginos are preferred from these considerations. 

^This upper limit on A is somewhat arbitrary but fairly unimportant in the present case since A grows very 
quickly when it becomes larger than unity. 
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Figure 2: Region in the rrih-Tn^ plane that is consistent with unification (shaded region) and 
precision electroweak data (below the dashed and sohd lines). 



In the model under study there are no dimension-five operators, so the dominant decay 
mode is expected to be p ^ vr^e"'". The high unification scale obtained above means that the 
proton will be long-lived regardless of the details of the unification mechanism. If proton decay 
proceeds via SU{5) gauge bosons, the decay rate is given by 

nP - vrV) = p^(l + D + Ffal [A^ + (1 + IK.HX] (18) 

where /tt = 0.131 GeV is the pion decay constant, D = 0.81 and F = 0.44 are chiral Lagrangian 
factors, Oat is a coefficient related to the n^p operator matrix element, and and are cor- 
rection factors due to the running of the couplings. A recent lattice-QCD calculation gives j2S] 
|aAr| = 0.015(1) GeV^, where the uncertainty is purely statistical. The systematic uncertainty 
is probably much larger; we take it to be ~ 50% [2^. The correction factors A^'^ split into 
long and short distance pieces: Ai^pt = Ai Y[^=i where Ai comes from the renormalization 
group evolution below Mz and Af'^ from that above. Here, Ai ~ 1.3 is identical to the SM 
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value, while the short distance pieces, to one-loop order, are 
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Using Mg = 2.8 x 10^^ GeV, and a^^ = 35.1, we find 

t{p 7r°e+) = 3 X lO^^^Vrs, (20) 
well in excess of the Super-Kamiokande bound of r(p — 7r*^e^) = 5.3 x 10^^ yrs |21j . 

3 Flavour Mixing 

Extending the matter content of the SM also introduces new sources of fiavour mixing. With 
mirror quarks, this pattern can be quite complicated, involving right-handed couplings to the 
W, and tree-level fiavour- changing couplings to the Z and the Higgs. We consider first the 
generic case, taking the most general set of Yukawa and mass terms possible. Next, we simplify 
our results making use of the fact that, in the model under study, the mirror quarks couple 
only to the third generation quarks and calculate explicitly the couplings of the heavy quarks 
to the weak gauge bosons and the Higgs. In subsequent sections we shall use these results to 
investigate the Higgs phenomenology, the collider searches for mirror quarks, and CP violation 
in B (f)Ks decays. 

Let A„ and be the fiavour-space mass matrices describing the fiavour mixing between 
gauge eigenstates. These matrices will be 4 x 4 and 5x5 respectively, and will have contributions 
from Yukawa couplings and Dirac mass terms. Both matrices can be diagonalized by bi-unitary 
t r ansf ormat ions : 

X^ = U^D^Wl Xa = UdD,Wl, (21) 

where the f/'s and W^s are unitary, and and are the diagonalized mass matrices. The 
corresponding (unprimed) mass eigenstates are then related to the (primed) gauge eigenstates 
by the unitary transformations 



AB„,B 



U 



L ; 
PQjQ 



R 
iP 
R 



AB„,B 



U 



(22) 
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Here, the indices A,B — 1,...4 correspond to {u,c,t,x}, while P,Q — 1,...5 refer to 
{d, s, b, cu, ^} respectively. 

In terms of the physical (mass eigenstate) fields, the charged currents become 



_ 7A» t 

'w- ~ '^W+ ' 
where the 4x5 flavour-mixing matrices are given by 



VS^ = Wrwf. (24) 



The matrix Vl is analogous to the CKM matrix Vckm of the SM. It is nearly unitary in the 
sense that VlVl^ = I^xa and V^Vl = Isxs — Vd, where the matrix Vd is defined below. The 
matrix Vr describes right-handed couplings, and has no analogue in the SM. 
Similarly, the hadronic neutral current is 

cos^^J^ = u^Y{\-l^^^^Oy,)ui + uiY{~sm^ey,)ui 

+dh'{~l + I sin' d^.)d^L + sin' 

+ I {uirvru'n - d'nrvf^di + F.^Vf^d^) (25) 

where the matrices Vu,Vd,Vd are given by 

(26) 





" u " u 











The off-diagonal elements of these matrices describe FCNCs. Each is Hermitian and satisfies 
3.1 Heavy Quark Neutral and Charged Currents 

The expressions above can be simplified considerably by using the fact that, in the model under 
study, the (gauge eigenstate) mirror quarks couple only to the quarks of the third generation (see 
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Eq. (PJ). The mixing between the mirror quarks and the first and second generation quarks is 
thus very small, and so will be neglected. Moreover, the mixing between the SM quarks is given 
approximately by the usual CKM description. The flavour violating effects among the heavy 
quarks are then related to the mixing of the right-handed SU (2) quark-doublet with the third 
generation right-handed quarks, as well as the mixing of the left-handed SU{2) quark-singlet 
with the left-handed bottom-quark. 

The mixing in the top sector must be very small to avoid a conflict with the 13{b 57) 
predictions j2ZI- We shall therefore assume, for simplicity, vanishing mixing in the top sector 
(1/4 = 0). The top-sector mass matrix is then diagonal, with rrit = -^Ut and = Mi. 

Following [Oj we take 1/5 = as well. The bottom-sector mass matrix, in the basis (6', cu', ^'), 
is then given by 

/ n F3 \ 

\d=\ Y2 Ml (27) 

V M2 / 

where Yi = -^Vi, i = b,2, 3. The phenomenologically interesting regime is Y^, Y3 <^ Y2, Mi, M2 [HI- 
Working to linear order in the small quantities Y^/Mi and Y3/M2, the left- and right-handed 
mixing matrices are 

' clCl Sl clSl 
Uh. = \ -sl cl I (28) 
-Sl Cl 



and 



W, 



Cr Sr 
-Sr Cr 
1 



(29) 



where sr = sm6R, sl = soiOl-, and sl = suiOl are given by 

Y2 



Sr 
Sl 
Sl 



Y,Y2 

Y^ 
M2 



Applying the mixing matrices to A^, the b-sector masses are 



(30) 



Y2 \ -1/2 



Mo 



(31) 
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To obtain a good agreement between the predictions of the model and precision electroweak 
data the angle in the right-handed sector must be sizeable, 

tan^i^ = Fs/Mi =^ 0.7, (32) 

while sl should be small, 

sin^L~ 0.09. (33) 

Note that Eq. ()32|) fixes sl in terms of the b mass. 

In this approximation, the relevant neutral currents read 

c| sin^ ew\ (T . , u \ ^R^R 



( Co sin (^w\ n n , \ 
^Rl^^R ( -y + 3 J - {hBTi^'^R + h.c) 



,„ , c| sin^6'vy 



2 3 



2 



1 sin^ 6'w\ ^ ,, . sin^ 6^^ 



+ -.7^-. + ^ j + ^uY^R-^- (34) 
Within the same approximation, the charged currents read 

= hYi^LbL + SLUJL + SL^L) 

+ XrI^ (crujr - SRbn) , (35) 

where we have neglected terms of order ml/Mf. 

In the above we have neglected the effects. We have verified that even large values of 
produce a negligible effect in the relevant right-handed Zbb coupling, while inducing only a very 
small change in the left-handed Zbb coupling (that can be compensated by a minor change of 
2/3), and so Y2/M1 = 0.7 still gives the bit fit to electroweak data. On the other hand, since y^ 
induces a mixing between the two heavy mirror quarks, the Z couplings and the mass spectrum 
of the mirror quarks are modified by this coupling. In particular, for non- vanishing values of 7/5 
there is always a mirror quark state in the b-sector with mass smaller than m° = (M^ + 1^2^)^/^. 
Therefore, for any given value of m^, mJJ, provides an upper bound on the lightest mass of the 
down-type mirror quarks. With Y2/M1 = 0.7 and Mi < 250 GeV as required by precision 
electroweak data the lightest down- type mirror quark can not be heavier than 300 GeV. 
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3.2 Higgs Couplings 

One of the most important immediate goals of high energy physics is to understand the mech- 
anism of electroweak symmetry breaking. In the Standard Model and its supcrsymmetric 
extensions, this symmetry is broken spontaneously through the vacuum expectation value of 
one or more scalar Higgs bosons. The same is true for the model under study and it is therefore 
quite relevant to understand the possible modification of the Higgs boson search channels at 
the Tevatron and the LHC. 

In addition to introducing new sources of flavour mixing, mirror quarks also modify the 
couplings to the Higgs. The Dirac mass terms for the mirror quarks mean that the Higgs- 
quark couplings need no longer be flavour diagonal in the basis of mass eigenstates, nor be 
proportional to the quark masses. We find that the coupling of the Higgs to the b quark is 
reduced relative to the SM. This, along with the contribution of the heavy quarks in loops, has 
interesting consequences for the detection of the Higgs. 

As in the previous section, we will assume that the mirror quarks mix almost exclusively 
with the third generation quarks and take ^4,^/5 ~ 0. This implies that the only Higgs-quark 
coupling that is significantly modified from the SM is that of the b-quark. The relevant terms 
in the Lagrangian are therefore 



Using the value tan^^ = Y2/M1 ~ 0.7 favoured by the model, we find that the hbb couphng is 
reduced by a factor of ~ 2/3. 
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4 Higgs Phenomenology 



4.1 Higgs Production and Decay 

This scenario modifies the phenomenology of the Higgs in two ways. First, by reducing the 
hbb couphng by a factor of c^, the partial width r{h bb) is attenuated by the square of 
this factor, cj^ ~ 4/9. Since this channel is dominant for Higgs masses below rrih — 130 GeV, 
the reduction of the partial width for this mode increases the branching fractions of the other 
modes in this range. Secondly, u quark loop effects increase the partial width r{h gg). This 
increases both the branching fraction of this mode, and the Higgs production cross-section by 
gluon fusion. 

Let us consider the effect of the u quark in a bit more detail. The presence of this quark 
in a loop connecting the Higgs to two gluons modifies the h gg partial width. Neglecting 
light quark contributions, and keeping only the effects of the dominant Yukawa coupling 7/2, 
the partial width becomes 



(40) 



where the function Fi/2{Tg) is given by j2Hl 

Fi/2 = -2rJl + (l-r,)/(r)], (41) 

with r„ = 4(^)2 and 



^ -i[ln(W7;_)-m]2; r<l; = 1 ± v^T^ ^^^^ 



r. 



The corresponding expression in the SM is obtained by setting sr = 0. Since the new term 
interferes constructively, the effect is to increase the decay width. While the h gg mode isn't 
directly observable at hadron colliders, this decay width is nevertheless important because it 
determines the cross-section for Higgs production by gluon fusion; a{gg ^ h) oc T{h ^ gg), 
up to soft gluonic effects. The 77 decay width is similarly modified by an u loop. In this 
case, the new contribution interferes destructively with the SM terms, the dominant parts of 
which come from W and Goldstone boson loops. However, the change in r(/i 77) is very 
small for any reasonable input parameter values. Figure El shows the dominant Higgs decay 
branching ratios in the model under study. Additional NLO corrections to the h ^ gg mode 
were included as well, following 
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Figure 3: Higgs Branching Ratios with Mirror Quarks, for = 250 GeV and 12/^^1 = 0.7. 

Figures |3 and El show the enhancement of a few Higgs discovery modes relative to the SM, 
which is mostly due to the increase in the gluon fusion cross-section. For low Higgs masses, 

~ 150 GeV, there is an additional enhancement of certain modes as a result of the decrease 
in the h -^bb branching ratio. It should be noted, however, that such low Higgs masses worsen 
the fit to the precision electroweak data in this model. 
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Figure 4: Enhancement of low mass Higgs 
production and detection modes for = 
250 GeV, Y2/M1 = 0.7. 
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Figure 5: Enhancement of intermediate mass 
Higgs modes for = 250 GeV, F2/M1 = 0.7. 
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If the Higgs mass exceeds 150 GeV the process h VV, where V is a real or virtual vector 
boson, becomes the primary Higgs discovery mode at both the Tevatron and the LHC [HIH \'M\. 
The inclusive modes are enhanced by the increase in gluon fusion, while modes in which the 
Higgs is produced by other means, such as vector boson fusion, are very slightly attenuated 
due to Higgs decays into mirror quarks. 

We have checked that including 7/5 has only a small effect on these results. Like the Zbb 
couplings, the hbb coupling is only slightly modified by y^. The additional mixing due to has 
a larger effect on the mirror quark couplings, increasing the h^C, coupling while decreasing that 
for hujuj. These changes precisely cancel each other in Eq. pn|l to the extent that -Fi/2(Tt^) = 
-^1/2(^5), which holds to within 7% for quark masses greater than the Higgs mass. 

4.2 Higgs Searches at the Tevatron and the LHC 

The enhancement of Higgs detection signals decreases the collider luminosity needed to find 
the Higgs. We have translated the above results into collider units using detector simulation 
results from |^ ESI Ell Tevatron, and jSHESI for the LHC. Figure El shows the minimum 

luminosity per detector (with CDF and D0 data combined) needed for a 3-cr signal at the 
Tevatron. Figure U\ shows the luminosity needed for a 5-0" discovery at the LHC. All channels 
displayed in this plot are for CMS alone except for the WW and ZZ modes, which are for 
ATLAS. These plots correspond to inclusive searches unless stated otherwise: VV —>■ h —>■ 
X denotes weak boson fusion channels, while tth X and W/Zh X denote associated 
production channels. For the inclusive channels, we have assumed that gluon fusion accounts 
for 80% of the total Higgs production. 

The model significantly improves the likelihood of observing the Higgs at the Tevatron for 
a Higgs mass between 120 and 180 GeV. Note that gg ^ h ^ t^t~ becomes the dominant 
discovery channel at the Tevatron collider in the low Higgs mass region. 

The predictions of the model for the LHC are less dramatic, although the improvement 
in the gg —>■ h '-fj process will make a light Higgs much easier to see. For intermediate 
Higgs masses, the inclusive h WW channel becomes competitive with the VV ^ h WW 
channel. For masses larger than the ones displayed in the Figure searches for the Higgs at 
the LHC can proceed via the golden mode h ZZ. 
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Figure 6: Minimum luminosity needed for a 3-a Higgs signal at the Tevatron for rriu, — 250 GeV, 
Y2/Mi = 0.7. 
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Figure 7: Minimum luminosity needed for a 5-a Higgs discovery at the LHC for m^^ = 250 GeV, 
Fs/Mi = 0.7. 
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5 Mirror Quark Collider Signals 



If the model is to improve the electroweak fits, the mirror quarks must not be too heavy. In 
particular, this requires < 250 GeV, which (using Eq. (HZj)) implies < 300 GeV as well. 
On the other hand, is largely unconstrained. These relatively low masses suggest that the 
Tevatron may be able to see mirror quarks by the end of Run II. Previous searches for exotic 
quarks have concentrated on a possible fourth generation b' quark. In the most recent of these, 
CDF has put a lower bound on the b' mass of rrib' > 199 GeV jTH], provided the branching ratio 
b' bZ is 100%. This bound is relevant to our model as well. 

At the Tevatron, mirror quarks are produced mostly by qq annihilation, with a smaller 
contribution from gluon fusion. Previous calculations of the top-quark pair-production cross- 
section apply to mirror quarks as well. These indicate aqq ~ 3.0—0.5 pb, for = 200—300 GeV 
at the centre of mass energy ^/s = 2.0 TeV [321 EZl- This is small, but comparable to the top 
production cross section in Run I, att = 6.1 ± 1.1 pb, where we have averaged the results of D0 
and CDF 

The up-type x quark is most strongly constrained in the model. It decays almost entirely 
by Xr ~^ bjiW due to a large tree- level right-handed W coupling, Eq. (j3^ . This will produce 
a signature very similar to that of the top quark. Indeed, top quark decays present a nearly 
irreducible background. Searching for the x therefore reduces largely to a counting experiment 
in which one compares the number of measured top events to the number expected. Searches 
at Run I of the Tevatron have already put interesting limits on since the top production 
cross-section measured there agrees well with SM predictions jSZ| (See, however, the comments 
in |SH].) 

In order to make a quantitative estimate of the present bound on a possible sequential 
top quark, we conservatively assume equal acceptance and detection rates for both x and top 
events. The fraction of x events in the top sample will be 40-10% for 200 < < 250 GeV. 
(In practice, the detection efficiency for x's will be slightly better since the b-tagging efficiency 
improves with jet Pt-) Comparing the Run I value for the top production cross-section with 
the theoretical prediction [HZI, we get that the cross-section for any new sequential top quark 
should be lower than 2.9 pb, at the 2-a level. This method, based just on counting, leads to a 
bound of about 

> 200 GeV. (43) 
This in turn implies > 250 GeV, well in excess of the b' bound from |19j . 
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At Run II, the goal is to determine the top production cross section to an accuracy of 
7-9% with a few fb~^ of data. Assuming that this goal is achieved, and considering that 
the uncertainty is mostly due to systematic effects (in particular, the proper determination of 
the b-tagging efficiency) and therefore weakly dependent on small luminosity variations, the 
Tevatron Run II will be sensitive enough to rule out the presence of a sequential top quark 
with a mass smaller than 230 GeV. This will imply an indirect bound on m^^ > 285 GeV. On 
the other hand, the Tevatron might see evidence of a x quark with mass smaller than 220 GeV 
at the 3-0" level. 

That the x ^ bW vertex is {V + A), Eq. makes the x somewhat easier to find. This 
is because the l^"'"'s emitted in xr ~^ bfiW~^ have positive or zero helicity, whereas those from 

— * h^W^ have negative or zero helicity ^. Leptons emitted by positive helicity W^s tend 
to be harder than those from longitudinal or negative helicity VT's. Thus, a slightly higher 
lepton Pt cut will increase the relative acceptance of x's, although the improvement will be 
small since the majority of W^'s emitted are longitudinal. CDF has looked for positive helicity 
PF's in top decays. They find a positive helicity fraction of = 0.11 ± 0.15 jlO], consistent 
with both the SM, and a x of mass above about 200-250 GeV, for which we predict a value of 

< 0.08-0.02. 

Run II at the Tevatron will also cover part of the mass range of the uj quark by direct searches 
for this particle. The strong ub mixing leads to tree-level bZu and bhu vertices, Eqs. ()34|) . ()39 |) . 
with the same 0{1) flavour-mixing factors. The dominant decay modes are thus ujr —>■ bnZ, 
and u ^ bh provided the Higgs isn't heavier than the u. Other modes are suppressed by loops, 
small flavour-mixing factors, and in the case oi u ^ phase space. Indeed, this decay is 
forbidden for almost all of the model parameter space consistent with precision electroweak 
data [Oj. If the Higgs is heavier than the u, the CDF bound applies directly to Couj — > bbZZ 
modes and constrains the u mass to be greater than 199 GeV. 

Things are more interesting if the Higgs is lighter than the uj. In this case, the ratio of the 
decay widths of the Higgs and Z modes is [Hj 

- bh) ^ {i-r^Y 

V{u^bZ) {l-rzY{l + 2rz) ^ ^ 

where = {mh/m^^Y , rz = {Mz/m^S)'^. Figure |H1 shows the branching ratios for these modes 

for a Higgs mass of 170 GeV. In their b' search, CDF looked for b'b' bbZZ events in which 

■^Correspondingly, the M^^'s emitted in the charge conjugate decays have negative or zero hehcity for a 
{V + A) vertex, and positive or zero hehcity for a. {V — A) vertex. We shall refer to both positive helicity 14^+ 's 
and negative hehcity W~'s as "positive helicity" W's, and so on. 
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one Z decayed into jets while the other decayed into a pair of high transverse momentum [Pt) 
leptons • They only accepted events in which the reconstructed mass of the lepton pair lay 
within the range 75-105 GeV and at least two jets were tagged as b's. For low Higgs masses, 
below 150 GeV, this search strategy is sensitive to both ujuj bZbZ and ujuj —>■ bZbh events, 
since in the latter, the Higgs decays predominantly into a bb pair which mimics the hadronic 
decay of a Z. This also lends itself to modifying the search strategy to include four b-tags [Hj. 
Even more search strategies become possible if the Higgs mass exceeds 150 GeV as favoured 
by the model. Such massive Higgs bosons decay mostly into WW pairs (see Figure El) so for 
instance, ujuj bZbh could be distinguished by looking for events with four jets, at least two 
of which are b-tagged, accompanied by a pair of high Pt leptons and large missing Et- 

In order to estimate the possibility of observing an uj quark at the Tevatron using the CDF 
b' search strategy, we shall assume that the Higgs mass is rrih = 170 GeV, and that this search 
strategy has no sensitivity to the uj ^ bh modes and a detection efficiency of 13% (as in Run I 
for large my). The most important background comes from Z^ events associated with hadronic 
jets. To reduce the background, one can impose a cut on the total transverse energy of the 
jets [in]. For > 250 GeV, a cut of ^ Et > 150 GeV will eliminate most of the background 
without reducing the signal in a significant way. The number of observable signal events scales 
with the luminosity and is approximately equal to 

Nu^c ^ 1.5-7.0 £[fb-i] , (45) 

for an uj mass varying from 300 GeV to 250 GeV. Due to the smallness of the background, a 
simple requirement for evidence of a signal is that at least five events be observed. Therefore 
the Tevatron Run II should cover the whole mass range of the model, m^^ < 300 GeV, if the 
luminosity is above 4 fb For a luminosity of 2 fb^^, a signal may be observed up to masses 
of about 280 GeV. 

The isosinglet ^ decays predominantly by biZ, t^W, and ^ ^ bh (see 

Eqs. fn^ . (jHK|) and (jHTHl ). All three of these go at tree level. Figure IHl shows the corresponding 
branching ratios for these decays. The ^ ^ bZ mode is dominant for < 400 GeV, although 
the ^ — i> tW quickly becomes important above the tW^-pro duct ion threshold 02] • For relatively 
light ^'s, of mass less than 250 GeV, the search strategies are similar to those for the uj. The 
CDF b' search is thus sensitive to a light ^. Their b' mass bound also applies in this case, 
although this limit may be weakened somewhat depending on the mass of the Higgs. 

In the preceding discussion we have taken 1/5 = 0. We have checked that including has 
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Figure 8: Branching ratios for decays of the uj 
quark with ruh = 170 GeV. 



Figure 9: Branching ratios for decays of the ^ 
quark with ruh = 170 GeV (and 2/5 = 0). 



only a small effect on the x ci^nd uj decay signatures. The most significant change for these 
decays is that the relationship between the masses of the x and u quarks, Eq. (fT7|). need no 
longer hold, as discussed at the end of section ITTl The branching fractions of the uj are changed 
only slightly. Keeping has a larger effect on ^ decays, and tends to decrease the branching 
fraction for ^ tW relative to the other modes. 



6 CP- violation in B (pKg Decays 

The value of the CP- violation parameter sin (2/?) measured in 5° — »• (pK^ decays appears to 
disagree with the value extracted from J/ ipKs decays: [IHl HH ES] 



sin(2/3) 



+0.735 ±0.054; 
-0.39 ±0.41: 



(46) 



This discrepancy is particularly interesting because the B — > (pKg mode is loop mediated, 
making it much more sensitive to new physics than the B J/ ipK^ mode, which goes at tree- 
level. We investigate whether this discrepancy can be explained by the FCNC's which arise in 
the model. 

It is not sin (2/3) that is measured directly, but rather the time-dependent CP asymmetry 
acp{t). For decays of the B^i= bd) meson into a CP eigenstate /, this is defined to be 

TiB%t)^f)-TiB%t)-.f) 



''CP 



T(B%t)^f) + T{B%t)^f) 
Cf cos(AMt) - Sf sin(AMt) 



(47) 
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where AM is the mass difference between the mass eigenstates, and Sf, Cf are given by 

Cf = SIgf, S, = (48) 
with ^ f given by (for the B'^ system) 

In the SM, the B (pKg amphtude has the form A = XtAo, where Aq is CP-invariant and 
A* = Vckm^c'km- The phase of A* is very small ^46 , so to a good approximation C^k = 0, 
and S(f)K = sin(2/5). This result applies to the B — > J/ipKg mode as well. For this reason, the 
values oi Sf measured in the (pKs and J j'^K^ modes are sometimes quoted as sin(2/3), as we 
have done in 

Physics beyond the SM can change the values of S^Ks and Sj/^Ks by adding additional 
CP-violating terms to the decay amplitudes. We have investigated whether the new tree-level 
FCNC Z-couphngs 

Jz ^ {sLl'VtbL - SRYVfbn) + {h.c.) (50) 

2cost^i„ V / 

can explain the apparent difference between S(j,K and Sj/^Ks- To simplify the analysis, we have 
assumed that the mixing of the mirror quarks with the SM quarks (other than the b) is very 
small. In particular, we have neglected all right-handed W couplings and all FCNC couplings 
other than those connecting the h and s quarks. \V'^^\ and IV^^^^I will also be treated as small 
parameters whose size we will bound below. These assumptions imply that the three generation 
CKM description of flavour mixing in the SM is correct up to small modifications, and that we 
can ignore the loop contributions of mirror quarks to the effective sh vertex. 

6.1 Constraints from Semi-Leptonic Decays 

The strongest constraints on these couplings come from semileptonic b ^ s modes 07j. At 
energies much below Mz, the SM decay amplitudes can be written as the matrix elements of 
an effective Hamiltonian; 



K// = K//(6 - 57) - ^A,(C97g9, + CfJIfQioJ, (51) 



where Qg^ = (sb) (^y- A) (H) v, s^ndQio^ = (s6)(y_A)(^0 a are four-quark operators. At = V^''^]^,jV^^j^j^, 
and I = fi,e. (See jlH] for a definition of 7ieff(b — > 57).) The FCNC couplings contribute to 
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Qgy and Qio^, and generate the new {V + A) operators Q'g^ = and Qiq^ = 

{sb)(V+A)ill)A- 

We neglect the contribution of mirror quarks to 7ie//(& s'j), since these only arise from 
loops, that become negligible for 7/4 = 0, contrary to the dominant tree-level effects included in 
our analysis. The effective Hamiltonian thus becomes 

^eff = n^ffib S7) - ^A, + Ciogio. + C^gW + CioQ'ioJ • (52) 

ysb 

In terms of rj := -j- and //':= — 4-^ the Wilson coefficients are now given by 



C9y = Cf*^ + (i-2sin2^Jr/, C^^ = (i - 2 sin^ 
r' — nSM r^i _ i„/ 



(53) 



Since \ — 2sin^^^u, ~ 0.05, to a good approximation we need only consider the shift in 
the Cioa coefficients. We have examined the effect of modifying the Cio^ coefficients on the 
branching ratios of the inclusive B XJ^l^ mode, as well as the two exclusive B —>■ Kl^l^ 
and B K*l^l^ decays. 

The constraint on b —>■ s FCNC's from the B —>■ Xsl^l^ mode has been considered previ- 
ously {e.g. IIHI im EOJ EI]). We repeat the analysis for this particular model using updated 
experimental results. From |1H]; modified to include the new (V^ + A) operators and neglecting 
lepton masses, the shift in the branching ratio relative to the SM is 



AB^^ = B-B^^' (54) 



a" 



cb 



{\CioJ' + \C[,J'-\C^oy) Bib^ceu), 



where Q := ^d, f{z) = 0.54 ± 0.04 is a phase space factor, k = 0.879 ± 0.002 is a QCD 

-ySM 
'10 a 



correction, and B{b ceu) = 0.109 ± 0.005. We have also taken Cf*^ = -^^^ ~ -4.2, 



a ^{nib) = 129, and 



Vts 



2 



1 in our analysis. The measured branching ratio and the SM 



prediction for this mode are listed in Tabled as is the 2-a allowed shift in the branching ratio 
based on these values. 

The inclusive modes B — > Kl^l^ and B K*l^l~ have been considered in j5()[ 15 Ij. The 
result for the latter mode, neglecting lepton masses, is [HH! 

AB^' = (4.1+^:?)xlO-«(|CioA-C^;oAp-|C^itT) (55) 



+ (0.9;°;^)x10-«(|(7ioa + C';oJ^-|C 



2 I r^SM 1 2 
10a I 
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Mode 






2-a Allowed Range (10~^) 


B Xsl^l' 




4.2 ± 0.7 54J 


-1.1 < AB^^ < 5.5 


B Kl+l- 


0.76+^:1^ [55, 56] 


0.35 ±0.12 [54^ 




-0.09 < AB^ < 0.91 


B K*l+l- 


1.68+^;^^ ± 0.28 56 


1.39 ±0.31 [54 




-1.3 < AB^' < 2.1 



Table 1: Experimental results and SM predictions. All errors were combined in quadrature, 
and the SM predictions were averaged over fi and e modes. 

Table HJlists the experimental results, the theoretical SM prediction, and the corresponding 2-a 
allowed range for AB^*. 

For the B Kl^l~ mode, the shift in the branching ratio is [20] 

where tb = 1.60 ±0.05 ps is the total B lifetime, and / is an integral of form factors. Explicitly, 
I = dsXfis)flis), where ^ ^ < S < and Ak(s) = 1 ± ± _ 2S - 2r^s 

with rx = (f^K/^B)'^ ■ We have evaluated the integral / numerically using the form factors 
in and find I = {0.056to^l^)- Again, table [T]hsts the relevant input data. 

Figure ^1 shows the 1] and rj' values consistent with all three semileptonic decay mode 
constraints taken at the 2-a level. Note that points within the two regions are correlated. 
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Figure 10: 2-a allowed ranges of 1] and rj' . 
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6.2 Range of S(jjK 

As for the semi-leptonic modes, the non-leptonic B (pK^ decay amphtude can be written in 
terms of an effective Hamiltonian. In the SM, this is given by [SZj 

10 

Xu{ClQl + C2Q^^) + Xc{C^Q'i + C2Q',)-XtJ2'^^Qi , (57) 



njSM 
^eff 



G p 

7i 



where A„ 



^CKM^CKM products of CKM factors, Qi{^j) are four-fermion operators, Cj(/i) 



are the Wilson coefficients, and /i is the renormahzation scale. ^ The operators Qi, Q2, Qs, • • • Qe-, 
and Qt . . . QiQ are the usual SM current-current, QCD penguin, and electroweak (EW) penguin 
operators respectively, as defined in jSIj. Note that the four-quark operators Qg and Qio are 
different from the semi-leptonic operators Qg^ and Qio^ considered above. 

If we include the FCNC couplings from the mirror quarks, the tree-level contribution to 
Ti-eff at scale Mw is 



AH, 



eff 



r]{sb) 



(V-A) 



(58) 



where the sum runs over q = u,d, s, c, b, g^^j^ = (T3 — sin^ 9^) is the Z{qq) l^r coupling, and 
7] = and rj' = —-^ are the same as above. The first operator, multiplied by r/, can be 
written as a linear combination of Qa, . . . , Qiq. The second operator, multiplied by 77', has 
no SM counterpart. We introduce a new "(l^ + A)" operator basis Q'^,...,Q'^q related to 
Q15 • • • , Qio by the interchange {V — A) {V + A) wherever these appear. 

We incorporate the new {V — A) operator contribution by modifying the Wilson coefficients 
at scale Mw- The changes are 



Cl^{Mw) 



Cl''{Mw) + 

2 

{Mw) + - sin^ e^T], 
Cl'\Mw)-\{l-siYiH^)'n. 



(59) 



For the {V + A) operators, [Cf^^)' {Mw) = 0, while the FCNC contribution gives 

1 , 



C',{M, 



w 



6 



V 



(60) 



^In writing "Heff in this form, we have made use of + Ac + Xt 
mirror quarks are included. 



1. This is also approximately true when 
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sin^ 6^7]', 



with all others zero. 

The RG evolution of these operators proceeds much like in the SM since the {V — A) and 
{y + A) operators evolve independently. The anomalous dimension matrix that determines 
running of both the {V — A) and the {V + A) operators is the same as in the SM. This follows 
from our definition of the (\^ + y4) operators, and the fact that these are renormalized by parity- 
invariant gauge interactions. We calculated the Wilson coefficients at scale = 2.5, 5.0 GeV at 
one- loop order in both the QCD and QED corrections using the results of [HHI- To this order, 
the corresponding initial values of the Wilson coefficients are taken at tree-level in the QCD 
corrections, although we have included the one-loop electroweak corrections which give a large 
contribution to C(){Mz). (This agrees with the conventions of |SZlEHlEn|-) 

Hadronic matrix elements for the B — > ipK^ transition at scale /i = 2.5, 5.0 GeV were 
estimated using factorization. Following jHOl IHI] , the amplitude is given by 



where Aq = —\/2Gpf(j,mfj,Ff^{m'^){e* ■ px) is a CP-invariant product of form factors and 
constants, and the Oj are functions of the Wilson coefficients. To leading order, they are 
o-2i-i = C'?i-i + Tr~C2i, and = €21 + tt— (^24-1, where N^ff is an effective number of 
colours. In writing (jUT|) . we have neglected annihilation contributions. While these may be 
significant, they do not introduce any new weak phases and are less than 25% of the SM penguin 
contribution to the amplitude jHIj. We find that the variation of the penguin coefficients with 
Neff is considerably larger than this. 

The Oj coefficients were calculated numerically by taking as input the 2-cr allowed values of 
T) and T)' from the previous section, and running the Wilson coefficients down to fi = 2.5, 5.0 
GeV. From these we calculate S^k, and the B —>■ (pKg branching ratio. This helps to reduce 
the theoretical uncertainty due to the sensitivity of the amplitude to variations of N^ff. Using 



(jnH) and the input parameters = 0.233 GeV, Ff^ = 0.39 ± 0.03, = 1019.4 MeV, 
rriK = 497.7 MeV, tub = 5279.3 MeV, tbo = 1.54 ±0.02 ps, |At| = 0.040 ±0.003, the branching 



A{B (f)K) = AoXt 03 ± 04 ± 05 - - {aj + ag + Oio) - -(a'^ ± ± a'^o) 



(61) 
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ratio is 

B{B^ (j)K^^) = (5.13 ± 0.15) X 10"^ 



as + a4 + as - ^{aj + ag + aio) - ^{a^ + a'g + a'^g^ 



(62) 

The range of S^k obtained for Nejf = 2, 3, ... 10 is shown in Figure [HJ and is plotted against 
the branching ratio. CLEO, BABAR, and BELLE have recently measured this branching 
ratio inaESEll, and the average of their results is i3(5° ^ 0ir°) = (7.7 ± 1.1) x 10"^ 
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14 16 
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Figure 11: Range of S^k accessible by the model plotted against the branching ratio for two 
choices of the renormalization scale. The vertical dotted hues indicate the 2-a allowed region 
for the branching ratio. 

We find that a range 0.2 < S^k ^ 1-0 can be explained by FCNC's in this model while 
simultaneously accommodating semi-leptonic B decay and B (pK branching fraction data at 
the 2-a level. While the shift in S^k from FCNC's is not large enough to completely explain 
the current experimental value, it is still significant, and reduces the discrepancy to below 2- 
a. A strong phase S in the new physics relative to the SM would only decrease the range of 
<S(j)K inn] • Setting 6 = n gives a result similar to that displayed in Fig. [TTl The result shows 
a strong dependence on the renormalization scale, fi, as well as the effective number colours. 
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Neff due to sensitive cancellations between terms in the amplitude. While this situation would 
be improved by adding higher order corrections, we do not expect such terms to change these 
general conclusions. 

The above result is more constraining than the one obtained in Ref. in which the 
effect of a vector-like pair of singlet down quarks on iS^x was considered. As we have done 
here, these authors investigate the range of S^k that can be obtained from the sZb vertex 
that is induced by vector-like down quarks. Our results should reduce to theirs in the limit 
rj' ^ 0, which corresponds to considering only the flavour mixing effects due to the singlets. 
By the same reasoning, the range of SfpK should be greater in the present model since even 
more flavour mixing is possible. Instead, these authors find a larger range for S^k than we 
have obtained. One of the possible sources of discrepancy between our analysis and the one 
in Ref. [^^1 resides in our use of a more stringent quantitative analysis of the constraints 
coming from the semileptonic B-decays. Another source appears to be their inclusion of the 
"colour-suppressed" operators (s/3&a)(y-A)(saS^)(y±^) at tree-level. While such operators do 
arise from QCD corrections to the sZb vertex, all tree-level effects may be described by the 
"colour allowed" operators {si3bi3)(v-A){saSa)(v±A)- 

7 Conclusion 

We have investigated the phenomenological properties of Beautiful mirrors, an extension of the 
Standard Model consisting of additional vector-like "mirror" quarks with the same quantum 
numbers as the SU{2) quark doublet and down quark singlet in the Standard Model. These 
exotic quarks mix with the bottom quark resulting in a modified value of the right-handed 
bottom quark coupling to the Z gauge boson, in agreement with indications coming from the 
precision electroweak data. A good fit to the precision electroweak data also demands that the 
additional quarks have masses lower than about 300 GeV implying a rich phenomenology at 
the Tevatron and LHC Colliders, as well as a possible impact on the CP-violating observables 
measured at the B-factories. In addition, the unification of gauge couplings is greatly improved 
within the model. 

In this article we have provided a detailed analysis of the question of gauge coupling unifi- 
cation. We find that the gauge couplings unify at Mq = (2.80 ±0.15) x 10^^ GeV. Perturbative 
consistency and stability of the model restrict the possible values of the masses of the Higgs 
and the mirror quarks. The allowed range, ruh = 170 ± 10 GeV, overlaps with the range of 
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values of these parameters which give the best fit to precision electroweak data [H] . 

Flavour mixing due to the mirror quarks leads to right-handed Z couplings, a very small loss 
of unitarity of the CKM matrix, and FCNCs. The flavour mixing also modifies the coupling 
of the b and uj quarks to the Higgs, while the couplings of the other quarks to the Higgs are 
not changed significantly. This has some interesting implications for Higgs searches at the 
Tevatron. In particular, the required luminosity for a Tevatron or LHC Higgs discovery in the 
WW decay mode, as well as in the t~^t~ mode at the Tevatron and the 77 mode at the LHC, is 
greatly reduced within this model. We have analyzed the search for mirror vector quarks at the 
Tevatron collider, and have found that Run II with a total integrated luminosity of about 4 fb~^ 
will be able to test all of the mirror quark mass range consistent with electroweak precision 
data. Finally, the b ^ s FCNCs which arise in this model can help explain the discrepancy 
between the values of sin(2/5) measured in the B (pK and B J/ipK decays. 
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Note added : As this work was being completed, Ref. [HZI appeared, in which the effect of 
tree- level {V ± A) sZb couplings on S^k was calculated. The modifications to the electroweak 
scale Wilson coefficients as well as the final result in this paper agree with our analysis. 
After this work appeared, revised experimental values of S^^k were reported by the Babar and 
Belle experiments. The new world average is still more than two standard deviations away from 
Sji^Ksi although there is also a. 2.1 a discrepancy between the two experimental results 
A related work, [^2] has also recently appeared, in which the search for Q = —1/3 vector-like 
quarks is investigated. The conclusions of these authors for the Tevatron are similar to ours. 
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